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Abstract. We investigate the extrinsic spin Hall effect in the electron gas model due 
to magnetic impurities, by focusing on Ce- and Yb-impurities. In the dilute limit, the 
skew scattering term dominates the side jump term. For Cc-impurities, the spin Hall 
angle asH due to skew scattering is given by — 87r<5 2 /7, where 6 2 (<SC 1) is the phase 
shift for d (I = 2) partial wave. Since asu reaches O(10 _1 ) if 82 £ 0.03, considerably 
large spin Hall effect is expected to emerge in metals with rare-earth impurities. The 
present study provides the highly efficient way to generate a spin current. 
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1. Introduction 

The spin Hall effect (SHE) has been attracting a great deal of interest due to its 
fundamental as well as technological importance. The intrinsic SHE originates from 
the Berry phase of the multiband Bloch function, so that the spin Hall conductivity 
(SHC) takes a material-specific value that is independent of impurity scattering in the 
low resistivity regime [TJ[2]. In transition metals, the intrinsic SHC takes a considerably 
large value 0, H], because of the "effective Aharonov-Bohm (AB) phase" induced by 
(i-orbital degrees of freedom and the atomic spin-orbit interaction (SOI) [7J [8]. 
Intrinsic SHE in Sr 2 Ru04 had been calculated in ref. [5] based on the t 2g tight-binding 
model. By the same mechanism, SHE in transition metals had been studied based on the 
realistic multiorbital tight-binding model in ref. [7J and the first-principles calculations 
[9]. The obtained SHCs agree well with recent experimental observations, in both the 
magnitude and sign [8]. 

In addition to the intrinsic mechanism, SHE is also induced by impurity scattering 
with the aid of SOI; this phenomenon is known as the extrinsic SHE [TOT [TTl [12]. It 
consists of a skew scattering term (cx s s H oc p _1 ) and a side jump term (<7g H oc p°); This 
fact suggests that considerable large SHC can be realized due to the skew scattering 
term in the very low resistivity metals, beyond the material-specific intrinsic SHC. For 
this reason, study of the extrinsic SHE is significant for realizing excellent spintronics 
device to produce or detect the spin current. 

The extrinsic SHE due to nonmagnetic impurities or vacancies had been studied in 
refs. [TOT [TTl fT2] . by developing the theory of the extrinsic anomalous Hall effect (AHE) 
[T3T [HI [15]. The obtained SHC is, however, very small unless the enhancement factor 
for the SOI due to the multiband effect is taken into consideration. For this reason, it is 
not easy to estimate the magnitudes of <7g s H and o"g H due to nonmagnetic impurities. In 
transition metals, there is no apparent experimental evidence for a large extrinsic SHE 
due to natural randomness. 

Dilute magnetic impurities such as Fe and Ce have been found to induce a large 
extrinsic AHE in simple metals [TTl EE] • The Hall angles in AuFe [T7] and rare-earth 
doped metals [16] reach 0.01~0.001 under high magnetic fields. Since the SHE has 
a similar mechanism to that of the AHE, one might expect that a giant SHE can be 
realized by magnetic impurities. However, theoretical analysis of the extrinsic SHE in a 
simple s-electron metal due to magnetic impurities has not been performed previously. 

In this paper, we study the extrinsic SHE based on a single-impurity Anderson 
model for Ce and Yb atoms, both of which are typical magnetic impurities. 
Nonperturbative effect of the large SOI for /-electrons (~ 3000 K) is correctly accounted 

f 2 |e| 

for. Due to skew scattering, the obtained spin Hall angle asH = — , which is the 

fy h 

ratio between the transverse spin Hall current and the longitudinal electric current, 
exceeds 0.1 for 5 2 > 0.03, where <5 2 (^C 1) is the phase shift for d (/ = 2) partial wave; 
It is more than 10 times greater than that in nonmagnetic metals such as Pt [I]. Since 



Giant Extrinsic Spin Hall Effect due to Rare-Earth Impurities 



3 



the spin Hall angle represents the efficiency of the creating spin current, giant spin Hall 
current is expected to emerge in simple metals by introducing dilute magnetic impurities. 
Moreover, the extrinsic SHC can take a large value well above the Kondo temperature 
Tk, in contrast to the anomalous Hall conductivity (AHC). The origin of the giant SHE 
is found to be the large SOI and angular momentum in the rare-earth impurities, which 
had not been treated appropriately in previous studies of the extrinsic SHE. 



2. Model and Hamiltonian 



In this study, we use the single-impurity Anderson model for Ce and Yb atoms. In a 
Ce 3+ ion with a 4/ 1 configuration, the J = 7/2 level is about 3000 K higher than the 
J = 5/2 level due to the strong atomic SOI [18]. Therefore, we consider only J = 5/2 
states for the Ce-impurity. In the same way, we consider only J = 7/2 states in Yb 3+ 
ion with 4/ 13 configuration. Here, we introduce the following single-impurity Anderson 
model for Ce/Yb atoms with both d- and /-orbitals [16J: 

H = E SfccLcfar + E E d dl m d am + E E f flf M 

k,a kam kM 

+ E {VLAadcm + h.C.} + E {VkMa^kafM + h.C.} 
kam kaM 

+ ~y E M M'i (!) 

Z M^M' 

where we have employed the electron (hole) picture for the Ce (Yb) impurity. c ha is 
the creation operator of a conduction electron with spin a = ±1. y M is the creation 
operator of a /-electron with total angular momentum J = 5/2 (7/2) and z-component 
M (—J < M < J) for Ce 3+ (Yb 3+ ). d\ m is the creation operator of a <i-electron with 
angular momentum m(— 2 < m < 2) [TjJ]. = k 2 /2m is the energy for the conduction- 
electrons, and (E d ) is the localized /(d)-level energy. Vj[ Ma . and V£ m are the mixing 
potentials, which are given by 

vL* = ^V f J2a% (T Y™(k), (2) 



V£ m =V4nV d Y 2 m (k), (3) 

where a^ a is the Clebsch-Gordan (C-G) coefficient and Y; m (A;) is the spherical harmonic 
function. We will show that the phase factor in Y; m (fc) and are indispensable to 
realize the SHE. For J = 5/2, = -a {(7/2 - Ma) /7} 1/2 8 m ,M-°/2 PHI EI]. Note 
that the SOI was neglected in the study of the AHE in ref. [16j. In the present study, 
we neglect the crystalline electric field of the /-orbitals, since it is small due to the small 
radius of the /-orbital wave function [16]. We put h = 1 hereafter. 

To discuss the scattering problem, it is useful to derive an effective Hamiltonian 
for the conduction electrons by integrating out the / and d electrons in eq. (CD). The 
obtained Hamiltonian is given by [IB] 

H c = E ^cLcfca + E J k,k'c{aCk'a + E J L,k'a' C ka C k'a' , (4) 



ka k,k'a kk'a,a' 
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where 



j d — v d V d * 

° k,k' rp d / j r km y k'm 

H — Hi m 



= 47rJ d J2Y 2 m (k)[Y 2 m (k')}\ (5) 

m 

jf - 1 V V s V s * 

J ka,k'a' ~ fit V kMo V k'Mo' 

H -El M 

= ^Jf E a™*%«>Y 3 m (k)[Yr'(k')]*, (6) 

Mmm' 

and J d = \V d \ 2 /(/j - E d ) and J f = \V f \ 2 /(fi - E f ). Here, & = E* + Re£ f ; is 
the /-electron self-energy due to the Coulomb interaction [20J. According to the 
scaling theory [22J, E* approaches the Fermi level as the temperature decreases due 
to the Kondo effect: \Jf\ is strongly enhanced near Tk, and below Tk, J/iV(0) ^> f 
due to strong resonant scattering, where N(0) = rnkp/^ir 2 is the density of state of 
the conduction band per spin. We also assume that J d N(0)(= — tan5 2 /7r) <C 1 since 
\H-E d \ ~0(leV) [16]. 

3. Calculations 

3.1. T -matrix due to c-f Resonant Scattering 

In this section, we study the scattering problem. In the both models given by eq. ([T]) 
and eq. (j4j), the T-matrix due to the c-f resonant scattering is equivalent; it is given by 

fel,<7l 

where its diagrammatic expression is shown in Fig. [I] (a), iV is the number of fc-points, 
and Gfc = (e — £fc) _1 - We note that the term containing both j{. ak i a i and J^ fc / given in 
Fig. [T] (b) vanishes identically due to the orthogonality of spherical harmonic functions. 
The solution of eq. ([7]) for k = k' and a = a' is simply given by 

where we have used the relations Y.mWL„? = 3\V f ? and i Ew V/ M „G° = 
|V/| 2 ^(£)5mm'- g{s) = jfJ2k^k( £ ) IS the local Green function. Assuming approximate 
particle-hole symmetry near /z, we put g R (0) = g(+i5) = — i7cN(0). Then, the 
quasiparticle damping rate in the T-matrix approximation is given by 

nN(0)I 
'1 + (irN(0)J f 

where n imp is the impurity concentration. Note that eq. is exact for n imp <C 1. 



7/ = -^im P ImT^ fc(T (0) = 3n imP i ; ^ ATl J T ^ , (9) 



3.2. Skew Scattering Term 

In this section, we study the skew scattering term using linear response theory. Initially, 
we consider the case N(0)Jf <C 1 and J d <C Jf, where the Born approximation is valid. 
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Figure 1. (a) Diagrammatic expression for the T-matrix due to c-f resonant 
scattering, (b) The diagram that vanishes after k\ -summation. 



In analogy to [T5], [16] , the skew scattering term is given by 
a SH - 2n n ^v N 2 2-, 2 dk x dk> 1 h{ jl 

k,k'a x V 

x l^(o)| 2 {r£ 2 £(o)4, fc + c .c}, (io) 

where — e (e > 0) is the electron charge. Its diagrammatic expression is shown in Fig. 
1(a). Here, 

k 1 ,a 1 

= 9 R (0)jL, k >*Jf (11) 

is the second order term of the T-matrix; the first order term in J^ tT ^ (r does not 
contribute to cr||j up to the first order term in J d [16]. Note that any diagram that 
contains the part shown in Fig. [1] (b) vanishes identically. The retarded Green function 
is given by (0) = (// — + ry) _1 , where 7 represents the quasiparticle damping rate. 
Here, we put 

7 = 7/ + 7o, (12) 

where 70 is the damping rate due to nonmagnetic scattering, such as c-d scattering 
(jd = 57mi mp iV(0) Jj) and the scattering due to disorders. The charge current is given by 
j° = —edek/dk v , where v = x, y. The spin current is then given by f v = (a/2)d£k/dk u . 
First, we consider the angular integration in eq. (fTUl) . which is given by 

^ 2 ( jka.k'gJk'.kfik, feQ) n ) (13) 
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Figure 2. Diagrammatic expressions for SHC induced by skew scattering (a) within 
the lowest order (extended Born approximation) contribution, and (b) in the T-matrix 
approximation with full order diagrams. 



where /(fe, k') 



m 2 ctefc dey 
k 2 dk x dk' y 



sin 9k cos (pk sin 6y sin <f>y , and (■ • -) n denotes the aver- 



age over the Fermi level, which is defined as [A(ki, ■ ■ -, k 



dOfcj • • • dVlk n 



• • fen). Since f(k, k') = 2vr {^(fe) - Yf^k)} [Yf^kf) + Y?(k')} /3i, angular integra- 
tion such as / dVtkY^ 1 a ^ 2 (k)Y l m (k)Y{ izl (k) appears in eq. ffTB"]) . This integral is finite 
only when I = 2,4. Therefore, the interference of the /(Z = 3) and d (/ = 2) partial 
waves is essential for skew scattering [16J. After performing the angular integrations, 
eq. (JTBl is given as 



eq.([T3D 



4 



Using the relations |Gf(0)| 2 » - e*) for small 7, and ^E^Gj^O) = g R {0) 

— Z7riV(0) , eq. fflOl) is transformed into 

= - ^—n imp J d J 2 f N(0) 
Since 7/ = 3nn imp N (0) J 2 in the Born approximation, af n is given by 



(14) 



a. 



SH 



27r2l7r2 J ^ 7 



.4 7/ 



(15) 
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(16) 



Now, we derive the skew scattering term using the T-matrix approximation, which 
gives the exact result for n imp <C 1. In this case, T^ 2 ^(0) in eq. ffTUj) is replaced with 

the full T-matrix T J ^ fe/o .(0) = T/^ 2 ^(0) (l — g R (0)Jf\ , where the first order term in 

Jkak'a nas been dropped. The diagrammatic expression for crg^ is shown in Fig. [2] (b). 
The angular integration in the T-matrix approximation can be performed as 

V2 \ ka > kaK k ' k 'I si 7l + ( 7 riV(0)J / ) 2 

Thus, cjgjj in the T-matirx approximation is given by (1 + (7rN(0)Jf) 2 ) 1 times eq. (JTj 
Therefore, the expression for the skew scattering term in eq. ([151) is valid beyond the 
Born approximation, by considering 7/ given by eq. (GO). 

Here, we discuss the temperature dependence of af K . From eq. (fT5l) . a^ K is 
proportional to 1/7/ when 7/ 3> 70- Since 7/ oc — In (T/Wband) for T ^> Tk [18J, 
°sh decreases as the temperature decreases. In contrast, af n cx 7^ for 7^ 70, which 
increases at low temperatures if 70 is constant. 



3.3. Side Jump Term 

In this section, we briefly discuss the side jump term <7g H . In the T-matrix 
approximation, Og J H is given by [15] 



a, 



sn 



^Jf2 2^ °^n——\G (0)1 



X 



fc(T, fc'cr' 

T/* G*(0) + <i?~A> 



(17) 



whose diagrammatic expression is shown in Fig. [3J In contrast to the skew scattering 



rpfR 

1 k'o'ka 




2 UxJ ka,k'a' 



+ c.c. 




u y J k'a',ka 



+ C.C. 



Figure 3. Diagrammatic expression for the side jump term in the T-matrix 
approximation. Both two diagrams give the same contribution. Note that Jd is not 
necessary for the side jump. 



term, is not necessary for the side jump term. As reported in ref. |21j . the large 
anomalous velocity v% Mrj , which is not perpendicular to the Fermi surface, arises from 
the ^-derivative of the phase factor in the c-f mixing potential as follows: 



dV kM* a 
dk T 3VkM ° 



a 

-i M-- 



yf 

2 J kl + kl kMa ' 



;is) 
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The final expression for the side jump term in the T-matrix approximation is given by 

^ = (19) 

where we have used the relations J2m<t Ma\V kM(T \ 2 = 3\Vj\ 2 (l — sin 2 ^) and 
Em(j °" 2 |VfcA/<rl 2 = 6|V/| 2 . By using eqs. ([15]) and (1TUI) . the ratio ag j H /cx§ s H is given 

by — . Therefore, will exceed a™ in dirty metals. However, the skew 

y 37iJ d n SH 

scattering contribution will be dominant in simple metals with magnetic impurities 
for n imp < 0.01 (p < 10/^cm) p]. 

4. Discussions 

4-1. Origin of the Skew Scattering Mechanism 



In section 13.21 we have studied the SHC due to skew scattering by using the Green 
function method. In this section, we discuss the origin of the skew scattering mechanism 
based on the Boltzmann transport theory. For this purpose, we study a simplified 
two-orbital model with M = ±5/2, assuming the strong crystalline electric field. 
In this model, V[ Ma = -a^A^V f {^Y^{k) 5 M ^ /2a + Y^(k)5 M>5/2a }. Since 
V kMa oc Y^ 3a (k) oc e~ 3la ^ k approximately, the second order term of is simply 
given as T^ 2 ^ R ~ — iN(0)e~ 3la ^ k ~^ k '\ In the Boltzmann transport theory, the spin 
Hall resistivity due to skew scattering is p s s H oc Y^a a \k x k'w{ka — > k'a)) , where 

w represents the scattering probability, which is proportional to n imp T ka k , a + J k k , 
due to Fermi's golden rule in the present model. According to ref. [16], skew 
scattering occurs when the scattering probability includes an assymetric component 
w ss (k — > k') oc Ime^^^'' oc (k x k') z . In the present model, w ss arises from 
the interference of the / and d scattering channel, w ss G (rl^l^J^y + c - c ) • ^ n 
fact, w ss (ka — > k'a) oc Ime WJ ^ fc_<?;>fe '* ) in this model since J kk , contains the term 
Y^ 2 (k) Y-2 ±2 (k') oc e ±2l ^ k ~^ k '\ In summary, a conduction electron with a hybridizes 
with l z = —3a state due to the strong SOI, and therefore the spin-dependent skew 
scattering probability w ss (ka — > k'a) oc Ime" 7 ^ -9 ^') arises from the interference of 
the / and d angular momenta. Thus, the origin of the SHE due to skew scattering 
mechanism is well understood based on the simplified two-orbital model. 

4-2. Estimations of the Spin Hall angle and the Skew Scattering Term 

of H 2e 

First, we estimate the magnitude of the spin Hall angle tanagn = z-. The 

2 

en 

longitudinal conductivity is given by a xx = , where n = kp/3n 2 is the density 

of the conduction electrons. Then, the spin Hall angle is given by 

tana H = — J d N(0)^. (20) 
7 7 
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Since tan5 2 = — KN(Q)Ja [18], eq. (1201) can be rewritten as tanasH = &2 if 

7 7 

|Jd|iV(0) <C 1. In the case where 7/ ^> 70, «sh reaches ~ 0.3 for 5 2 ~ 0.1 |24j. 
Therefore, giant SHE will be realized in metals with rare-earth impurities. 

Recently, giant SHE was observed in FePt/Au devices [23]. It was reported that 
the large spin Hall angle in Au is ~ 0.1. Since the intrinsic SHE in Au is small due 
to the Fermi surface being mainly composed of s-electrons [7], the observed giant SHE 
may be explained by the skew scattering mechanism. Very recently, the extrinsic SHE 
in AuFe metal has been studied with the first-principles band structure calculations 
independently of our study [25] . 

Second, we estimate the magnitude of the skew scattering term for 8% = 0.1, where 
the spin Hall angle is 0.3. Then, <jg s H = (ft/2e)0.3 x o xx from the above discussion. In 
the case where the resistivity is 10/ificm, we obtain af u = 15000(ft/e)n _1 cm _1 if we 
put the length of unit cell a = 4 A, which is the case for Pt and Au. The obtained skew 
scattering term is about 60 times larger than the SHC observed in Pt: 240 (ft/ e)Q~ 1 cm~ 1 . 
Since the resistivity in Pt and Au is p ~ 10/ificm and p ~ lpQcm, respectively, the 
skew scattering term due to magnetic impurities should dominate the intrinsic term. In 
contrast, the skew scattering due to "natural randomness" has not been observed in the 
AHE in transition metals, even in p ~ 10/ificm |26j. Therefore, the orbital degrees of 
freedom of the impurity is indispensable for the giant extrinsic SHE. 

Finally, we compare the magnitude of af n and a S g H . If we put kp = n/a and 

7 « 7/ in eq. ([19]), of H = (4/3) x lO^/e^^cm" 1 for a = AA. Therefore, 

the skew scattering term will be comparable in magnitude to the side jump term when 
p ~ lOOyuficm. Thus, the skew scattering term is dominant over the side jump term 
when p <C 100/ificm. 

5. Summary 

In summary, we have studied the extrinsic SHE due to magnetic impurities, where Ce 
and Yb impurities have been discussed as a typical case. The obtained expression for 
the skew scattering term shows that the spin Hall angle reaches O(10 _1 ), which is more 
than 10 times larger than that in Pt. Therefore, we propose that the extrinsic SHE due 
to rare-earth impurities is useful as a method for creating large spin current efficiently 
The present study also suggests that the giant SHE will arise from other rare-earth 
atoms as well as Ad and 5d atoms that possesses strong SOI. As in the case for the 
intrinsic SHE, orbital degrees of freedom are essential for generating a huge extrinsic 
SHE. According to eq. (j20~|) . giant spin Hall angle always appears when 7^ ^> j Q . Thus, 
the Kondo resonance, N(0)Jf 3> 1, is not a necessary condition for the giant SHE. In 
contrast, the AHC (= a xy /H z ) is enhanced by the Kondo effect since it is proportional 
to the uniform susceptibility. 

In a similar way to the present study, we can obtain af K and Og J H for J = 7/2, 
which is the case for Yb-impurities. In this case, a^ a = {(7/2 + Ma) /7} 1 ^ 2 5 nij M-a/2, 
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and 7/ is given by (4/3) times eq. The obtained af H and <7g H are given by (—3/4) 
times eq. f fT5|) and eq. (fl~9l) . respectively. 

Acknowledgements 

We are grateful to D. S. Hirashima and J. Goryo for fruitful discussions. 
References 

[1] J. Sinova, D. Culcer, Q. Niu, N. A. Sinitsyn, T. Jungwirth, and A. H. MacDonald, Phys. Rev. 
Lett. 92 (2004) 126603. 
S. Murakami, N. Nagaosa and S.C. Zhang, Phys. Rev. B 69 (2004) 235206. 

E. Saitoh, M. Ueda, H. Miyajima and G. Tatara, Appl. Phys. Lett. 88 (2006) 182509. 
T. Kimura, Y. Otani, T. Sato, S. Takahashi, and S. Maekawa, Phys. Rev. Lett. 98 (2007) 156601; 

L. Vila , T. Kimura, and Y. Otani, Phys. Rev. Lett. 99, 226604 (2007); Y. Otani et al., 
(unpublished) 

H. Kontani, T. Tanaka, D. S. Hirashima, K. Yamada, and J. Inoue, Phys. Rev. Lett. 100, 096601 
(2008). 

H. Kontani, M. Naito, D. S. Hirashima, K. Yamada, and J. Inoue, J. Phys. Soc. Jpn. 76 103702 
(2007). 

T. Tanaka, H. Kontani, M. Naito, T. Naito, D. S. Hirashima, K. Yamada, and J. Inoue, Phys. 

Rev. B 77, 165117 (2008). 
H. Kontani, T. Tanaka, D. S. Hirashima, K. Yamada, and J. Inoue, arXiv:0806.0210. 

G. Y. Guo, S. Murakami, T.-W. Chen, and N. Nagaosa, Phys. Rev. Lett. 100, 096401 (2008). 
M. I. Dyakonov and V. I. Perel, JETP Lett. 13 467 (1971). 
J. E. Hirsch, Phys. Rev. Lett. 83, 1834 (1999). 

S. Takahashi and S. Maekawa, J. Phys. Soc. Jpn. 77, 031009 (2008). 
J. Smit, Physica 24 (1958) 39. 
L. Berger, Phys. Rev. B 2 (1970) 4559. 
A. Crepieux and P. Bruno, Phys. Rev. B 64 014416 (2001). 

A. Fert, J. Phys. F 3, 2126 (1973), A. Fert and A. Friederich, Phys. Rev. B 13, 397 (1976). 
A. Fert and O. Jaoul, Phys. Rev. Lett. 28, 303 (1972). 

A. C. Hewson, The Kondo Problem to Heavy Fermions (Cambrige Studies in Magnetism, 1993); 

K. Yosida, Theory of Magnetism (Springer 1996). 
We omit the atomic SOI for d-electrons since it is much smaller than \/i — E \. 
K. Yamada and K. Yosida: Prog. Thcor. Phys. 53 (1975) 970. 

H. Kontani and K. Yamada, J. Phys. Soc. Jpn. 63, 2627 (1994). 

F. D. M. Haldane, Phys. Rev. Lett. 40, 416 (1978). 

T. Seki, Y. Hasegawa, S. Mitani, S. Takahashi, H. Imamura, S. Maekawa, J. Nitta, and K. 

Takanashi, Nature 7, 125 (2008). 
Recently first-principles band structure calculation has been performed by the present authors, 
and §2 is estimated as ~ 0.1 for Yb impurities in Ag host. 

G. Y. Guo, S. Maekawa, and N. Nagaosa. larXiv:0809.29"2ll 

T. Miyasato, N. Abe, T. Fujii, A. Asamitsu, S. Onoda, Y. Onose, N. Nagaosa, and Y. Tokura, 
Phys. Rev. Lett. 99 (2007) 086602. 



